Wavepacket dynamics in energy space, RMT and quantum-classical correspondence 
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We apply random-matrix-theory (RMT) to the analysis of evolution of wavepackets in energy 
space. We study the crossover from ballistic behavior to saturation, the possibility of having an 
intermediate diffusive behavior, and the feasibility of strong localization effect. Both theoretical 
considerations and numerical results are presented. Using quant al- classical correspondence (QCC) 
considerations we question the validity of the emerging dynamical picture. In particular we claim 
that the appearance of the intermediate diffusive behavior is possibly an artifact of the RMT strategy. 



We are interested in the dynamics that is generated 
by a Hamiltonian of the type H — E + 'W where E is a 
diagonal matrix whose elements are the ordered energies 
{En}, with mean level spacing A, and W is a banded 
matrix. It is assumed that W is similar to a 'Banded 
Random Matrix' (BRM), with non- vanishing couplings 
within the band < |n — m| < b. These coupling ele- 
ments are zero on the average, and they are characterized 
by the variance a = ((|W„mp))^^^. Thus, there are four 
parameters (A, 6, cr, h) that controls the dynamics. One 
important application of BRM is in solid-state physics for 
the study of localization in quasi-one-dimensional disor- 
dered systems. In this frame non-zero values of A reflect 
the presence of a constant electric field along the sam- 
ple. However, in this Letter we mainly have in mind the 
original motivation following Wigner Namely, the 
study of either chaotic or complex conservative quantum 
systems that are encountered in nuclear physics as well 
as in atomic and molecular physics. For this reason the 
above defined model (with non-zero A) is known in the 
literature as Wigner's BRM (WBRM) model. 

Consider a system whose total Hamiltonian is Ti-iQ, P), 
where (Q, P) is a set of canonical coordinates. We assume 
that the preparation and the representation of the sys- 
tem are determined by a Hamiltonian 'Hq{Q, P). We also 
assume that both HoiQ^P) and Ti{Q,P) generate clas- 
sically chaotic dynamics of similar nature^. We choose 
a basis such that the quantized Hamiltonian matrix Tio 
has a diagonal structure TLq = E. According to gen- 
eral semiclassical arguments 0, the quantized Hamilto- 
nian matrix 7i, in the same basis, has a band structure 

= E -t- W. The WBRM model can be regarded as a 
simplified local description of the true Hamiltonian ma- 
trix. However, there is one feature that distinguishes the 
effective WBRM model from the true Hamiltonian. It is 
the assumption that the off-diagonal elements are uncor- 
related, as if they were independent random numbers. In 
this Letter we would like to explore the consequences of 
this RMT assumption on the dynamics. Below we define 
the classical limit of the WBRM-model, and the various 
parametric regimes in the quantum-mechanical theory. 
We analyze the dynamical scenario in each regime, and 
we explain that the emerging picture is incompatible with 
the quantal-classical correspondence (QCC) principle. 



Taking 7i(Q,P) to be a generator for the (classical) 
dynamics, the energy Tioit) — Ho{Q{t), P{t)) fluctuates. 
The fluctuations are characterized by a correlation time 
T^i, and by an amplitude 6E^,. The three parameters 
(A, b, a) that deflne the effective WBRM model are de- 
termined by semiclassical considerations H) . One obtains 
A oc W^, and b oc and a oc fi^'^~^^^'^, where d is 

the number of degrees of freedom (dimensionality) of the 
system. In this Letter, we flnd it convenient to define^ the 
two classical quantities (r^,, 5E^i) in terms of the common 
quantum-mechanical parameters: 



SE,, = 2Vba 



(1) 



The classical dynamical scenario is formulated by using 
a phase-space picture Q. The initial preparation is as- 
sumed to be a microcanonical distribution that is sup- 
ported by one of the energy-surfaces of HoiQ^P)- For 
t > 0, the phase-space distribution spreads away from 
the initial surface. 'Points' of the evolving distribution 
move upon the energy-surfaces of 7i((5,P). We are in- 
terested in the distribution of the energies 7io(0 of the 
evolving 'points'. It is easily argued that for short times 
this distribution evolves in a ballistic fashion. Then, for 
t 3> T^i, due to ergodicity, a 'steady-state distribution' 
appears, where the evolving 'points' occupies an 'energy 
shell' in phase-space. The thickness of this energy shell 
[Q equals 5E^i. Thus we have a crossover from ballis- 
tic energy spreading to saturation. The dynamics in the 
classical limit is fully characterized by the two classical 
parameters t^, and SE^^. 

We are going to study the corresponding quantum- 
mechanical scenario. We want to explore the conse- 
quences of assuming that WBRM model can be used as 
an effective model for the true dynamics. At t = the 
system is prepared in an eigenstate of Hq — E. For t > 
the evolution of the system is determined by 7Y = E-FW. 
The evolving state is ip{t), and we are interested in the 
evolving distribution |(n|-0(t))p. We shall use the fol- 
lowing terminology: The standard perturbative regime 
is (cr/A) < 1; The Wigner regime is 1 < (cr/A) < 6'^^; 
The ergodic regime is <C (c/^ ^ b^'^', The localiza- 
tion regime is <C (ct/A). See It is easily verifled 
that the limit h —i- corresponds to the ergodic regime 
or possibly (provided d = 2) to the localization regime. 
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The structure of the eigenstates aofH has been stud- 
ied in We denote the average shape of an eigen- 
state as Wsir) — {\(pa{na + r)]"^) where ipa{n) — {n\a), 
and Ua is the 'site' around which the eigenstate is lo- 
cated. The average is taken over ah the eigenstates that 
have roughly the same energy Ea ^ E. There are two 
important energy scales: One is the classical width of the 
energy shell 5E^^ , and the other is the range of the interac- 
tion Afc = bA. In the standard perturbative regime W-E{r) 
contains mainly one level, and there are perturbative tails 
that extend over the range Af,. In the Wigner regime, 
many levels are mixed: the main (non-perturbative) com- 
ponent of WE{r) has width F = 27r(cr/A)^ x A, and the 
shape within the bandwidth Ab is of Lorentzian type. 
However in actual physical applications this Lorentzian 
is a special case of core-tail structure where the tail 
can be found via first order perturbation theory. Outside 
the bandwidth the tails decay faster than exponentially 

. On approach to the ergodic regime We (r) spills over 
the range At,. Deep in the ergodic regime it occupies 
'ergodically' the whole energy shell whose width is SE^i. 
In actual physical applications the exact shape is deter- 
mined by simple classical considerations 1^,^. Deep in 
the localization regime W^ir) is no longer ergodic: A typ- 
ical eigenstate is exponentially localized within an energy 
range SE^ = £,A much smaller than SE^^. The localiza- 
tion length is ^ « 5^. In actual physical applications it is 
not clear whether there is such type of localization. To 
avoid confusion, we are going to use the term 'localiza- 
tion' only in the sense of having 6E^ <C SE^i. 

Now we would like to explore the various dynamical 
scenarios that can be generated by the Schrodinger equa- 
tion for an{t) — (ri|'0(t)). Namely, 

— - --E a --VW a (2) 

rn 

starting with an initial preparation a„ = Snm at t—0. 
In a previous study only the localization regime has 
been considered. Here we are going to consider the gen- 
eral case (A 7^ 0). We describe the energy spreading 
profile for /: > by the transition probability kernel 
Pt{n\m) = (|a„(i)p). The angular brackets stand for 
averaging over realizations of the Hamiltonian. In partic- 
ular, it is convenient to characterize the energy spreading 
profile by the variance M(t) = ^^^(n— m)^Pt(ri|m), and 
by the participation ratio N{t) = (X]n(-Pt('^l"^))^)~^' 
and by the total transition probability p{t), and by the 
out-of-band transition probability q{t). Both p{t) and 
q{t) are defined as X)! -^("-1'^) where the prime indi- 
cate exclusions of the term n = m or exclusion of the 
terms |n— to| < b respectively. Equation (j^) has been in- 
tegrated numerically using the self-expanding algorithm 
of to eliminate finite-size effects. Namely, additional 
106 sites are added to each edge whenever the probabil- 
ity of finding the 'particle' at the edge sites exceeded 
10^^^. Fig.l illustrates the time-evolution of the en- 
ergy spreading profile. From such plots we can define 



various time scales. The times tbaii and i^^t pertains to 
M (t) and mark the departure-time from ballistic behav- 
ior and the crossover-time to saturation. The time t^t^ 
pertains to N{t) and marks the crossover to a station- 
ary distribution. The time scale t-p^t pertains to p(t) 
and marks the disappearance of the simple perturba- 
tive structure (See below). The asymptotic value of 
q(t), if it is much less than 1, indicates that the system 
is either in the standard perturbative regime or in the 
Wigner regime, where out-of-band transitions can be ne- 
glected. The saturation profile is given by the expression 
P^{n\m) — '^a\{n\a)\'^\{a\m)\'^ , and it is roughly ap- 
proximated by the auto- convolution of Wsir). Therefore 
the saturation profiles (Fig.l) are similar to the average 
shape of the eigenstates. We have found that M(oo) sat- 
isfies a scaling relation Y = 2X-{1— exp(— 1/(2X)) where 
X=(cr/A)/6"/" and Y={M{oo)y^yb^. This scaling rela- 
tion is similar to the one that pertains to the average 
shape of the eigenstates [Q . 




FIG. 1. Left figure: Representative examples for the time- 
evolution of the energy spreading profile. The variance M{t), 
the participation ratio N{t), the total transition probability p{t), 
and the out-of-band transition probability q{t) are plotted for 
A=0.05,0.10, 0.20, 0.50 and 6=80. The units of energy are chosen 
such that (T=l and the units of time such that h=l. Right figures: 
Saturation profiles for various regimes: (a) Standard perturbative 
regime; (b) Wigner regime; (c) Ergodic regime; (d) Localization 
regime. The distance between the tick-marks on the horizontal axis 
of (a)-{c) is b. In (d) the full scale is |n-m|<2000. The full scale of 
the vertical log-axis is — 15< ln(P)<0. Note that the n=m term is 
factor 3 larger compared with its immediate vicinity [Q. In (a) and 
(b) the !/(?!— m)^ behavior of the (in-band) tail is fitted by dashed 
lines. Notice the appearance of a core region in (b), indicated by 
the 'flattening' of the profile for |n— m|<20. The in-band profile in 
(b) corresponds to a Lorentzian with a very high accuracy. 

In the standard perturbative regime each eigenstate of 
Ti is localized 'perturbatively' in one energy level. Thus, 
for arbitrarily long times the probability is concentrated 
mainly in the initial level. We can write schematically |^ 

Pt{n\m) sa (5„m + Tail(n-TO;t) (3) 

where 1&i\{n—m;t) — {a/h)'^tFt{{En—Em)/Ji) within 
the range of first order transitions (0 < |rt— m| < b), 
and zero otherwise. Here Ft{uj) — t-{smc{ujt/2))'^ is the 
spectral-content of a constant perturbation of duration t, 
and sinc(a;) — sin(a;)/x. We have trivial recurrences from 
n to m once t becomes larger than 2Trh/{En—Em.)- The 
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global crossover to quasi-periodic behavior is marked by 
the Heisenberg time in = 271%//^.. The total normaliza- 
tion of the tail is much less than unity at any time. 

In the Wigner regime, one observes that the pertur- 
bative expression (||) is still valid for sufficiently short 
times t ^ tprt. Let us estimate the perturbative break- 
time tprt- For short times {t < t^,) the spectral function 
Ft{aj) is very wide compared with the bandwidth Ah of 
first-order transitions. Consequently we can use the re- 
placement Ft {u) t and we get that the total transition 
probability is p{t) w & x [at/TiY . On the other hand, for 
t > r^i, the spectral function -Ft(a;) is narrow compared 
with the bandwidth, and it can be approximated by a 
delta-function. As a result we get p{t) w a'^/{hA) x t. 
The condition p{t) ~ 1 determines t^^^ leading to: 

_ J fiA/a^ for 1 < a/ A < Vb 

~ [ h/{Vba) for Vb < a/ A ^ > 

It should be noticed that for cr ~ A we get t^^^ t^- 
Thus, taking recurrences into account, we come again to 
the conclusion, that for cr <C A there is no perturba- 
tive breaktime. The variance {6E{t)Y = A^ x M{t) of 
the energy distribution (||) is easily calculated. We get a 
ballistic-like behavior, followed by saturation, 

5Eit)^\^^''t'j^^^' f'- (5) 
^ ' \ SE^i for t>T,i ^ ^ 

For t ~ iprt the tail (|^) becomes Lorentzian-like, and it is 
characterized by a width h/t = T. For t > tp^ expression 
losses its validity, but it is obvious that the energy 
cannot spread any more, since it had already acquired 
the saturation profile. 

It should be realized that neither (||), nor the 
Lorentzian-like saturation profile of the Wigner regime, 
could correspond to the classical spreading profile. In 
the latter case the saturation profile is characterized by 
two genuine quantum mechanical scales (F, A^), whereas 
the classical ergodic distribution is characterized by the 
single energy scale SE^i. See Fig. 2. However, in spite of 
this lack of correspondence, the variance ^ behaves in 
a classical-like fashion. Using the terminology of Q we 
have here restricted rather than detailed quantal-classical 
correspondence (QCC): The quantal Pt{n\m) is definitely 
different from its classical analog, but the variance SE{t), 
unlike the higher moments of the distribution, turns out 
to be the same. 

In the ergodic regime the time scale becomes larger 
than tprt, and therefore r^i loses its significance. At 
t ~ tprt the quantal energy-spreading just "fills" the en- 
ergy range A;,, and we get 5E{t) w A;,. The perturbative 
result (|) is no longer applicable for t > tp,t. However, 
the simplest heuristic picture turns out to be correct. 
Namely, once the mechanism for ballistic-like spreading 
disappears a stochastic-like behavior takes its place. The 
stochastic energy spreading is similar to a random-walk 
process where the step size is of the order A;,, with tran- 
sient time equals tp,t. Therefore we have a diffusive be- 
havior 6E{t)'^ — D^t where 



De = C-A^/tp,, = C ■ A'^b'^'cr/n cx n (6) 

and the numer ical prefactor § is C « 0.85. This diffu- 
sion is not of classical nature. The diffusion can go on as 
long as {D^ty^ < SE^i, hence the ergodic time is 

= b-''^ ha/A^ cx l/h (7) 

After this ergodic time the energy spreading profile sat- 
urates to a classical- like steady state distribution 
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FIG. 2. The upper diagram illustrates the various dynamical sce- 
narios which are described in the text. The flow of time is in the 
vertical direction. See for closely related diagrams. The lower 
plot illustrates the various energy-scales that characterize the as- 
sociated stationary distributions: The bandwidth A;, is indicated 
by a horizontal solid line; The width of the non-perturbative com- 
ponent is indicated by the grey filling; The width of the energy 
shell is indicated by the dashed line; The variance 5E(oo) is in- 
dicated by the bold solid line. In the localization regime we have 
SE{oo) RJ SE^ < SE^i. 

In the localization regime the quasi periodic nature of 
the dynamics is important. The 'operative' eigenstates 
are defined as those having a non-negligible overlap with 
the initial state m. These eigenstates are located within 
the energy shell whose width is SE^, . If the eigenstates are 
ergodic, then all of them are 'operative', and therefore the 
effective level spacing between them is simply A„ff « A. 
However, if the eigenstates are localized, then only ^ out 
of them have a significant overlap with the initial state m, 
and therefore the effective level spacing is A„ff sa SE^J^. 
The effective energy spacing A^ft is the relevant energy 
scale for determination of the crossover to quasiperiodic 
behavior. The associated time scale is t* — 2TTh/A^ti, 
and it may be either less than or equal to the Heisenberg 
time til = 21111/ A. The localization regime is defined by 
the condition t* < i„g- In this regime the diffusion stops 
before an ergodic distribution arises, and we should get 
D^t* « (5i?|. Inserting the definition of t* and solving 
for we obtain the well known [HJ] estimate fa b^. For 
the breaktime we obtain 



3 



t* = W'^n/a cx {i/h) 



(8) 



Note that the locahzation range is dE(^ 
If the diffusion were of classical nature, we would get 
SE^ oc [l/TiY^^ as in the semiclassical analysis of 
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FIG. 3. (a) The times tbaii, tsat, tata and tp^t are numerically 
determined. Different values of b are distinguished by the rela- 
tive size of the symbols. The axes are X = log((T/A)/ log(b), and 
Y = log{t/tH)/log{27r6). Note that Y=-l implies t=Td, and Y=0 
implies t=tH- In the ergodic regime t^at departs from tbai and an in- 
termediate diffusive stage appears. The saturation time approaches 
tn but eventually drops down once we enter into the localization 
regime, (b) Both log{M{oo))/ log(6) and log(A'"(oo))/ log(fe) arc 
plotted versus log((j/A)/ log(fe). The arrow indicates a global hor- 
izontal shift of the N{od) plot for presentation purpose. 

The various dynamical scenarios discussed above are 
summarized by the diagram of Fig. 2, and can be com- 
pared with the data presented in Fig. 3. As expected 
from the theoretical considerations we have in the Wigner 
regime U^u « t,at ~ t^i and i,ta ~ iprt- In the ergodic 
regime we have as expected itaii ~ iprt ^ '''ci, while 
^sat ~ ista = iorg- Thus, In the ergodic regime we have 
a premature departure of the ballistic behavior, and the 
appearance of an intermediate diffusive stage. 

Our major motivation for studying WBRM model 
comes from 'quantum chaos' (see introduction). Namely, 
WBRM model can be regarded as an effective model 
for the analysis of the dynamics of a 'quantized' clas- 
sically chaotic system. The condition to be in the regime 
(ct/A) ^ 6^'^ can be cast into the form h ^ Cp^t, where 
Cprt = SE^iT^i is a classical scale. In this regime the 
perturbative result Eq.(|^) is valid. The derivation of (|^) 
is not sensitive to the presence or the absence of subtle 
correlations between matrix elements. Therefore (^) is 
valid in case of the 'quantized' Hamiltonian, as well as 
in case of the effective WBRM model. Hence we may 
say that the applicability of an effective RMT approach 
is trivial in the regime h ^ Cp^t- In contrast to that, 
in the non-perturbative regime (h C^^t), correlations 
between matrix elements become important, and it may 
have implications on the dynamical behavior. Whether 



an effective RMT approach is valid becomes a non-trivial 
question in the non-perturbative regime. 

In the regime h ^ Cp,t we have restricted QCC [1). 
It means that QCC holds only for the variance SE{t). 
Fixing all the classical parameters, including the time t 
which is assumed to be of the order of r^i, we can al- 
ways define a sufficient condition h <C Csc for having 
detailed QCC Detailed QCC means that the quantal 
energy spreading profile Pt(n\m) can be approximated 
by a classical calculation. The considerations that lead 
to the determination of the classical scale Csc are dis- 
cussed in |5|. We cannot give an explicit expression for 
Csc because it is a non-universal (system-specific) pa- 
rameter. Detailed QCC implies that Eq.(|^) should hold 
again once the condition h <^ C'sc is satisfied. 

For the WBRM model we have found that for H <C Cp^t 
there is a pre-mature departure from ballistic behav- 
ior, followed by an intermediate diffusive behavior. So 
we have a contradiction here between RMT considera- 
tions on one hand, and QCC considerations on the other. 
Thus, if the RMT approach is non-trivially valid, then it 
is only in a restricted range Csc «C ?i <C Cp^f Out- 
side this regime it is either trivially valid {U ^ Cp^t) 
and we have restricted QCC, or else it is not valid at all 
{h <C Csc) and instead we have detailed QCC. It may 
be true that in many cases the RMT considerations are 
not valid for the purpose of analyzing time-dependent dy- 
namical scenarios. A similar situation may arise in the 
theory of quantum dissipation: There is one-to-one cor- 
respondence between the regimes in Fig. 5 of (|] and the 
regimes that have been discussed in this Letter. 
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Physically, going from Ho to TL may signify a change of an 
external field, or switching on a perturbation, or sudden- 
change of effective-interaction (as in molecular dynamics) . 
The numerical prefactors are chosen such that whenever 
^ applies, we have tbai ~ Td and SE{cx}) — SE^i. 
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